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The transition from a liquid to a glass in colloidal sus- 
pensions of particles interacting through a hard core plus an 
attractive square-well potential is studied within the mode- 
coupling-theory framework. When the width of the attrac- 
tive potential is much shorter than the hard-core diameter, a 
reentrant behavior of the liquid-glass line, and a glass-glass- 
transition line are found in the temperature-density plane 
of the model. For small well-width values, the glass-glass- 
transition line terminates in a third order bifurcation point, 
i.e. in a (cusp) singularity. On increasing the square-well 
width, the glass-glass line disappears, giving rise to a fourth 
order A4, (swallow-tail) singularity at a critical well width. 
Close to the A^ and A4, singularities the decay of the density 
correlators shows stretching of huge dynamical windows, in 
particular logarithmic time dependence. 

PACS numbers: 61.20.Ne, 64.70.Pf, 82.70.Dd 



I. INTRODUCTION 

Colloidal suspensions have been studied extensively be- 
cause of their practical importance and because of their 
relevance in biophysics. These systems are also of great 
theoretical interest since they are models for conventional 
matter. They can be prepared for a large span of den- 
sities so that the states can be gases, gels, liquids, crys- 
talline solids, or glasses. Light scattering can be used 
to measure the static structure factor and various cor- 
relation functions. The dynamics can be explored over 
a wide range of length scales and over huge dynamical 
windows . Fascinating with colloidal systems is that 
the interaction can be tuned to some extent by varying 
the coating of the particles and the composition of the 
solvent [Q-D. It is possible to realize the hard-sphere 
system (HSS), the basic model underlying all theories of 
simple liquids UJ^. One can also prepare systems where 
the hard core is complemented by an attractive shell. 
This allows to study the interplay of repulsion and at- 
traction. As a contribution to such studies, a theory for 
the glass formation resulting from a strong short-range 
attraction among densely packed hard-sphere colloidal 
particles shall be presented in this paper. 

In hard-sphere colloidal dispersions, the liquid-glass 
transition has been studied by van Megen and Pusey 



g. They measured correlation functions <j)q{t) for den- 
sity fluctuations of a representative set of wave numbers q 
over about four decades in time t. It was found that these 
correlations decay to zero as expected for a liquid only 
for packing fractions (p below a critical value (pc- At </3c, 
the long-time limit of the correlators, fq = (j)q{t — > 00), 
changes discontinuously to a certain value > 0, in- 
creasing further with packing fraction, fq is the Debye- 
Waller factor of the amorphous solid, i.e. of the glass, 
and generalizes the order parameter introduced by Ed- 
wards and Anderson in the theory of spin glasses §. The 
evolution of the glassy dynamics for the HSS was stud- 
ied comprehensively by van Megen and coworkers [p|-p^ . 
The data suggest, that it is the well known cage effect 
§ which causes the glassy dynamics and the arrest of 
density fluctuations at ipc- 

The cage effect is the essential physical concept under- 
lying the mode-coupling theory (MCT) for the evolution 
of glassy dynamics in simple systems fljjl^ . This the- 
ory allows the calculation of </'g(t) and thus fq from the 
equilibrium structure factor Sq. As a function of con- 
trol parameters like singularities of fq, called glass- 
transition singularities, may occur. The simplest type, 
called a fold bifurcation, describes a liquid-glass transi- 
tion at = It implies a subtle dynamical scenario, 
giving rise to universal features of glassy dynamics, which 
have been identified in a leading-order-asymptotic expan- 
sion of the MCT equations. A review of the basic results 
IS given m Ref. @. In Refs. detailed quantita- 

tive comparisons of the data for hard-sphere colloids with 
the MCT predictions are presented. It is shown that the 
theory accounts for the experimental facts within a 15%- 
accuracy level. An illuminating summary of these studies 
is given in Ref. Results for the shear modulus have 
also been interpreted with the universal MCT formulas 
JlSf . The evolution of glassy dynamics for ip increasing 
towards (£c was also studied for polymer micronetwork 
colloids |19[-p^. Here, the interparticle interaction is not 
known. But the authors demonstrated, that a consistent 
fit of their data with the universal MCT formulas was 
possible. Preliminary studies of the glassy dynamics of 
charge-stabilized colloids indicate, that these data can 
also be explained within MCT The reported find- 

ings shall be taken as a justification to base the theory 
in this paper on the MCT for simple systems. 



1 



Our studies deal with the square- well system (SWS), 
characterized by a hard-core repulsion for interparticle 
distances r < d, and by a constant attraction poten- 
tial within the shell o? < r < d -I- A. The theory fo- 
cuses on the high-density regime, say (p > 0.4, so that 
the cage effect is essential for the dynamics. The rel- 
ative attraction-shell width 6 — A/d is assumed to be 
small, say S < 0.15. The main outcome of our theory is 
the prediction of a higher-order glass-transition singular- 
ity at a critical packing fraction ip* somewhat above the 
critical point ipc of the hard-sphere system and a criti- 
cal width S* of about 0.04. This singularity organizes a 
subtle phase diagram and opens up various possibilities 
for glassy relaxation. The new results reflect the inter- 
play of two mechanisms for particle localization, i.e. for 
the arrest of density fluctuations. It can either be dom- 
inated by repulsion of the particle by its cage-forming 
neighbors, or by the formation of bonds to the bound- 
aries of the cage. Preliminary calculations ^4 25 1 based 
upon Baxter's adhesive-hard-sphere model 1 2^ hinted at 
some findings to be derived in this paper. Baxter's model 
treats the limit (5 — > 0, so it cannot deal with the indi- 
cated singularity at 6*. Moreover, taken literally, the 
Baxter model cannot be used as a basis of MCT appli- 
cations, since there appears a divergency due to excita- 
tions with large wave vectors. The results for this model 
in Refs. p^-[2q| depend in an ill-defined manner on the 
large-g cutoff used there, a problem which is avoided with 
the SWS. 

Dense systems of colloidal particles characterized by a 
hard core and strong attractions of a range smaller than 
the core diameter by a factor of at least 10 were realized 
experimentally, when adding nonadsorbing polymers to 
either a suspension of colloidal hard spheres [g8[ or to 
emulsions |^^, in solutions of sterically stabilized par- 
ticles when decreasing the solvent quality pl|-^, and 
in copolymer micellar systems when changing the tem- 
perature |3^. Such systems were also studied in Monte 
Carlo simulations 1 36 3^ . Non-equilibrium phenomena 
characterized by a number of aspects were found which 
cannot be understood from the glassy states formed in 
hard-sphere solutions. First, amorphous soHds could be 
formed by increasing the attraction strength even though 
the packing fraction was kept fixed well below the value of 
the hard sphere glass transition pO|-^. Second, increas- 
ing the strength of a short-ranged attraction by adding 
small polymers, melting of the glass states was reported 
for the colloid-polymer mixtures [^9112^ ] . Third, the non- 
decaying frozen structures seen when immersing polymer 
coated colloidal particles into solvents of decreasing qual- 
ity 1 32 1 exhibited a much larger Debye- Waller factor at 
small wave vectors than hard-sphere systems. This indi- 
cates a much higher rigidity of the solid states on inter- 
mediate length scales. In support of this observation, 
viscoelastic measurements for intermediate frequencies 
found strongly concentration dependent elastic moduli 
pO| , pT 33 3^. It will be shown that our results provide a 
qualitative explanation of the reported findings. 



The paper is organized as follows. In Sec. ||, we report 
our results for the structure factor of the SWS and discuss 
those features which cause various qualitative results of 
the MCT solutions. Section |l| presents the main result 
of this paper, showing the phase diagram, and discussing 
the properties of the glass states resulting from the in- 
terplay between attraction and repulsion. In Sec. IV, we 



present some results for the dynamics which illustrate 
that the higher order glass-transition singularities cause 
relaxation stretching which is much more pronounced 
than is known for the HSS. The last Sec.^ presents some 
concluding remarks. 



II. STRUCTURE FACTOR CALCULATIONS 

A. The Model 

The structure factor Sq is the essential input informa- 
tion needed to formulate the MCT equations. In this 
section, Sq shall be discussed for the square-well system 
(SWS). Only such states shall be considered for which 
Sq depends smoothly on the particle density p, on the 
temperature T, and on the wave number q. The inter- 
action potential V(r) for particles with separation dis- 
tance r consists of a hard-core repulsion for r < d, and 
it has the negative value —Uq within the attraction shell 
d < r < d + A. The structure can be specified by three 
control parameters: the packing fraction tp of the hard 
cores, the ratio 6 of thermal and attractive energy, and 
the relative width 5 of the attraction shell: 



ksT /uq, 6 — A/d . 



(1) 



Let us note the standard concepts needed for a discus- 
sion oi Sq 1^. g{r) and h{r) = g{r)~l abbreviate the pair 
distribution and the total correlation function, respec- 
tively. The Fourier transform hq of the latter determines 
the structure factor Sq — 1 + phq. The Ornstein-Zernike 
equation formulates an integral equation for h[r), where 
the kernel is the direct correlation function c(r). In the 
wave vector domain, it reads Sq = 1/[1 — pcq], where 



47r /■°° 
Cq ~ — / dr s\niqr)[rc{r) 

q Jo 



(2) 



Baxter's method of the Wiener-Hopf factorization [[4[|3q] 
shall be used to reformulate the Ornstein-Zernike equa- 
tion. The basic concept of this theory is the factor func- 
tion Q{r). It is defined as a continuous real function for 
r > 0, determining Sq via its Fourier transform: 



Sq' = Q{q)Q{qr , 
/>oo 

Q{q) = 1 — 2ttp / dr exp{iqr)Q{r) 
Jq 



(3a) 
(3b) 



It is anticipated that Q(r) as well as c(r) vanishes beyond 
a certain distance R. For < r < i?, there holds 
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rc(r) = -Q'ir) + 27rp / ds Q'{s)Q{s - r) . (4) 



Furthermore, one finds for r > 

rh{r) = -Q'{r) + 2Trp / ds{r ~ s)h{\r - s\)Q{s) . (5) 
Jo 

For the SWS, g{r) =0 is fulfilled for < r < d, and 
therefore, using h{r) — g{r) — 1, Eq. (||) splits into three 
subequations. Most simple is the result for the middle 
part, A < r < d, where the formula known from the 
theory for the hard-sphere system (HSS) is reproduced: 

Q'{r)^ar + h. (6a) 

Here, coefficients a and h are introduced by 

a = 1 — 27r/9 / dsQ{s), b = 2TTp / dssQ{s). 
Jo Jo 

(6b) 

Writing G{r) = rg{r), one finds for small distances, < 
r < A, 



Q'{r) =ar + b- 2Tip 



d+A 



dsG{s-r)Q{s), (6c) 



r+d 



and for the attraction shell, d < r < d + A, one obtains 



/•r — d 

Q'ir) ^ar + b- G{r) + 2ttp / ds G{r - s)Q(s) . 

(6d) 

Some approximation for c(r) has to be introduced into 
Eq. (^) in order to close the system of Eqs. (^ and (|^). 
In this paper, the Percus-Yevick approximation (PYA) 
and the mean-spherical approximation (MSA) shall be 
applied Q . Nezbeda already studied the structure factor 
for the SWS using the PYA for small well widths ^!^. 
His equations could be solved only in a restricted region 
of parameters. Since the boundary of this region of ap- 
plicability is close to the parameter region ip « 0.5, 9^1 
of interest in this paper, it does not seem appropriate to 
base the following calculations on these results. 



B. Approximations 

Within the PYA, one writes c(r) = g{r)[l — 
exp{V{r)/kBT)] outside the hard core. Substitution of 
this ansatz into Eq. (^) and using Eq. (6d) leads to the 
approximation for d < r < d + A 



g-«o/fcBTg(^) = ar + - 27rp / ds Q'{s)Q{s - r) 



d+A 



/>r — d 

+2ttp / dsG{r ~ s)Q{s) . (7) 
Jo 



Equations (^) and (|^) for Q(r) and G{r) are solved 
numerically. To proceed, the equations are discretized 
straightforwardly. On each of the three r-intervals, a 
grid of equally spaced points r„ is chosen, where n = 
1, 2, . . . , 1000. The functions Q'{r) and G{r) are calcu- 
lated iterating Eqs. (^ and (0). At each step, the func- 
tion Q{r) is evaluated from Q'{r) using a 5-point nu- 
merical integration. The procedure is carried out until 
the difference between two successive iterations summed 
over all points of the r-grid becomes less than 10~^^. The 
integral in Eq. ( ^ ) is determined by a simplified Filon 
procedure to obtain Q{q) and hence Sq. 

The MSA uses c(r) = -V{r)/kBT outside the hard 
core. Substituting this ansatz into Eq. (jj), after integra- 
tion one obtains for d < r < d + A 

rd+A 

Q{r) = 2TTp dsQ{s)Q{s~r) 

J r 

+ [{d + Af-r^]/i29). (8) 

Equations (||) and (^) are solved analytically in a lead- 
ing and next-to-leading order expansion, using the well 
width 6 as the small parameter. For the organization 
of the expansion, the quantity K — 5/9 is considered 
fixed. This procedure is motivated by Baxter's theory 
of sticky hard spheres |2^. He evaluated Sq in the limit 
5 — !■ 0, uq — > oo, keeping a parameter equivalent to K 
fixed. Details of the calculation can be found in the ap- 
pendix. The hard core diameter d shall be used as the 
unit of length. For ^ < r < 1, the factor function is the 
parabola known from the theory of the HSS: 



Qir) ^ ar^ /2 + br + c, 



(9a) 



with a, 6, and c now being smooth functions of the SWS 
control parameters. For < r < J, there is an enhance- 
ment above this parabola: 



Q{r) = ar^/2 + br + c + 2ipK'^ • 5 [1 - [r / 5)f 



(9b) 



Within the attraction shell, the leading order result de- 
scribes a linear decrease of Q{r) from K to zero. The 
leading correction adds a quadratic modification. One 
finds for 1 < r < 1 + (5 



Q{r) = K 




r-1 



1 - 



Here and in the following we denote the constants from 
Eq. (^ as a = Oo -I- K6 ■ ai, b — bo + K5 ■ bi, and c — 
cq + KS-ci. The leading order contributions are the result 
of the Baxter limit 5^0: 



ao 



l + 2ip 



(1 - ^)2 



12K(p 
(1^ 



(10a) 



(9c) 
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bo 



Co 



-3<y9 



2(1-^)2 

-1 



6Kip 



2(1-^) 



K . 



(10b) 
(10c) 



The terms in brackets exhibit the resuhs for the HSS 
Q. The coefficients of the next-to- leading-order contri- 
butions are 

ai = [6ip{5^ - 2) - 72co</?2(i _ ^)] /(^ „ ^ (^^^^ 

6i - [9(^(1 - 2ip) + 36co(p2(l - ifi)] /(I - ^)2 , (lib) 
ci = [1 - 7^ + 12co^(l - (^)] /(2(1 - ifi)) . (11c) 



Substitution of Eq. (0) into Eq. (3b) yields Q{q) as a 



combination of trigonometric functions. It is elementary 
to work out the somewhat lengthy expression and thus 
via Eq. ( |3a| ) the desired result for Sq. 

The large-g asymptote of the direct correlation func- 
tion Cq shall be obtained from Eq. (||) by evaluating the 
asymptote of the Fourier-sine transform of the function 
/(r) = rc{r). From Eq. (^ one concludes, that /(r) is 
smooth except for at most three points r^^^ = S, r^^-* = 1, 
and A^'' — 1 + 6. At these points there can be a dis- 
continuity, given by that of the derivative of the fac- 
tor function: = Q'(r(') - 0) - Q'(r(') + 0). Let 
us also note from Eq. (^) the initial value /(r = 0) = 
A = -Q'(0) - 6(/?(3(0)^ For the exact solution of the 
problem, A = must hold, but due to the approxima- 
tion scheme used here, a finite value of 0{S) remains: 
A = -K5 [12^co(ci + 2Kip) + bi] + ^(J^). Thereby, the 
Baxter resuh A = 0(5°), is improved. The can 
be determined easily from Eq. (^), in particular /*^^^ = 0. 
Integrating by parts, the integral in Eq. (||) becomes: 
[/^"^ + E» /^'^ cos((7r(*))] Iq + 0{l/q^). Hence one ar- 
rives at Cq = + 0{\/ q^), where the asymptotic tail 
reads 

c^^ = {An/q^) ■ [A + Bcos{q) 

+ 2C sin {q{l + 5/2)) sin {q5/2)] . (12) 

The second term in the bracket has a form familiar from 
the PYA result for the HSS. But the coefficient B is a 
smooth function of K and 5 which reduces to the HSS 
value for K = 5 = Q: 



B = a + b + K{l2ipco ~ 1) . 



(13a) 



The third term in the bracket is due to the existence of 
the attraction well. Its prefactor reads 



c = ii + 5)/e. 



(13b) 



C. Results 



The spinodal lines of the SWS are shown in Fig. |l| for 
three representative values of the well width 6. They 



specify the divergence points of the compressibility, i.e. 
the zeros of for q = 0. The spinodal is the bound- 
ary of the regime of absolute instability with respect 
to the liquid-vapor transition. Only states outside this 
regime can be considered in the following. Substitution 
of Eq. (H) into Eq. (|3b|) yields elementary expressions 
for 0(0) within MSA. We have not been able to deter- 
mine the spinodals within the PYA, due to numerical 
instability of the algorithm. The high-density regime in- 
vestigated in the following applications of the MOT is 
indicated as the strip between the two dotted vertical 
lines. 

Figure || exhibits structure factors Sq calculated within 
the MSA for S — 0.05, and the corresponding pair dis- 
tributions g{r), calculated numerically from Eq. (H), for 
states marked by diamonds in Fig. |l|. The S'^-versus-q 
curves exhibit a principal refraction peak as known from 
other simple liquids Q|. It is caused by the hard-core 
driven excluded volume phenomenon, g{r < d) = 0. The 
high-temperature curves 1 and 2 exhibit peaks, which 
are only slightly smaller and somewhat broader than the 
peaks of a HSS at the same densities. The attraction 
modifies the pair correlations and thus the excluded vol- 
ume effects, as can be inferred by comparing the curves 
1 and 3. Lowering T, the short-ranged attraction causes 
the particles to move closer, i.e. the most probable inter- 
particle spacing decreases. Therefore, the peak position 
shifts to higher q upon cooling. The distribution g{r) 
develops a more rapidly varying structure at distances 
which are multiples of the particle diameter, and this ex- 
plains the decrease of the peak height and the increase 
of the peak wings in Sq. A change of the density at low 
temperature modifies the peak in a similar manner as 
discussed above for large T, cf. curves 3 and 4. However, 
lowering (p drives the system closer to the spinodal, and 
therefore the exhibited change of Sq for small q is larger 
than expected for a HSS. 

Results for the pair correlation g{r) obtained by dif- 
ferent closures of Eqs. (^ and (|^) and by other meth- 
ods have been published by Lang et al. |^^. For both a 
small and a large well width considered there, S — 0.03 
and 6 = 0.5 respectively, we find our results for g{r) in 
agreement with the Monte-Carlo simulation results ob- 
tained by Lang et al. Only for the small well width, g{r) 
is underestimated in the well regime, 1 < r < 1 + S, 
by about the same amount as Nezbeda's approximation 
p9| , p0| overestimates g{r). At r > 1 + (5, however, our 
solution appears to be in better agreement with the sim- 
ulation results. This behaviour is similar to what holds 
for the optimized random-phase approximation pTf . 

The largc-g tail of Cq will be of importance in the fol- 
lowing. In Fig. ||it is shown that the asymptote, Eq. (p^), 
describes Cq very well for q > 20. The results have been 
evaluated for the state discussed in Fig. |^ with the label 3, 
where {A,B,C) = (-0.092,1.63,7). The taU consists 
of a part due to the first two terms in Eq. (p^, which 
differs from the HSS result merely by modifications of 
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the coefficients A and B. Tfiis part of the asymptote is 
shown in Fig. |^ (dotted line) in order to emphasize that 
the last contribution in Eq. can be dominant. The 
next-to-leading order contributions to our results are not 
relevant for a discussion of the qualitative features of the 
tail. Therefore, let us write the lowest order formula for 
the tail as c^^ = + c^**. Here, the coefficients of d^"^ 
are obtained via a Baxter-like limit, 5^0, K= 6/6 fixed. 
Noting A — > in this limit, we find 



^*usA = 0.5258, ei^sA = 0.2332, 



<5itisA = 0.0465. 



(16b) 



(47r/g2)Bo cos(g) 



(14a) 



It will be shown in Sec. Ill, that the two states specified 
above are of central importance. The values found for the 
corresponding densities (p* and well widths S* are close to 
each other. The difference in the values for the effective 
attraction strength 1/6* is well explained by Eq. (15). If 
one inserts « 7.5 as obtained from the MSA, one finds 
0.2233. 

All calculations within the MSA are based on the 



^MSA 



The attraction-induced tail in this approximation reads ^^^ll--^ expansion for the factor function, say^ Q^^^ir) 



c— = {An/q^) ■ K{2/5) sin(q5/2) sin(g(l + 5/2)) . (14b) 

For the example in Fig.|, [Bq, K) = (4/3, 1/3). For q be- 
low an upper cutoff q^ — t: /5^ the function in the bracket 
of Eq. ( |l4b[ ) increases almost linearly with q. This ex- 
plains the increasing importance of the c^*' contribution 
relative to the c^'^p one, as is demonstrated in Fig. ^. If 
K is sufficiently large or if |i?o| is sufficiently small, one 
can identify a lower cutoff qi = \Bq\/K such that c^'* 
dominates c^°P: 



d"^ = {ATlK/q) sm{q{l + 5/2)) , < g < q„ 



(14c) 



In the wave- vector interval between qi and g^, Cq ex- 
hibits a power-law decrease slower than the one of the 
true large-q tail, which dominates only for q^ q^- 

The PYA and the MSA differ solely by their ansatz for 
the direct correlation function c(r) within the attraction 
shell 1 < 7- < 1 + (5. Within the PYA, c(r) depends on 
r via the r-dependence of the pair distribution function, 
cpYA(r) = [l-exp(-uo/(fcsT))]g(r), while the MSA 
assumes a constant cmsa — ut)/{kBT). In this paper, 
systems with narrow attraction shells are of interest, 5 <C 
1. Therefore it is a reasonable approximation to ignore 
the r-dependence of cpya by writing g{r) ~ gd — g{r — 
1). Thus, for every state where a solution of the PYA 
exists, there is a solution of the MSA, yielding the same 
structure factor. However, the corresponding solution 
for the MSA has to be evaluated for an effective reduced 
temperature 6^^^ . The latter is a smooth function of 
6 = uo/kBT, ip and 5, estimated by 



1 



^MSA 



1 — cxp 



-Up 

knT 



9d- 



(15) 



Consequently, the PYA and the MSA yield the same sce- 
narios for the structure factor in the parameter regime 
of interest in this paper. This finding is demonstrated 
in Fig. ^ for the basic quantity of the structure factor 
theory, the factor function Q{r). The result calculated 
within the PYA for the parameter triple 



V?PYA = 0.5293, 6pYA = 1.1000, 



(5pYA = 0.0429 



(16a) 



is very close to the one obtained within the MSA for 



defined by Eqs. (||) to (y|. To control this result, 
Eq. (|) and the analog of Eq. (g) for the MSA clo- 
sure have also been solved numerically to get the cor- 
rect MSA factor function, say Q^^^(r). The difference 
5Q{r) = Q^^^{r) - Q°''P(r) is positive and about 3% 
(1%) for 5 = 0.25, (0.15) for r < 0.5; and it decreases for 
r increasing above 0.5. 

The dotted line in Fig. ^ exhibits the parabola for 
Q(r), Eq. (|9a|), for coefficients of the HSS. Introduction 
of the attraction smoothly renormalizes the coefficients 
a, b, and c, such that the parabola shifts upwards and be- 
comes flatter. There appears a region of positive values 
for Q{r) near the core surface r « 1. These shifts cause 
the smooth drifts of the Fourier transform for Q{q), which 
lead to the drifts of Sq discussed in connection with Fig. ^ 
and to the appearance of a spinodal, cf. Fig. |^. The only 
qualitative new feature, which is caused by the attrac- 
tion well, is the almost straight decrease of Q{r) within 
the interval 1 < ?■ < 1 + 5. Equation (^) yields the 
slope in leading order as Q'{r) = ~K/5. In the Bax- 
ter limit, this slope diverges. The specified almost con- 
stant part oiQ'{r) causes the attraction-tail contribution 
to t he as ymptote of c<j, Eq. (14b). The power law tail, 
Eq. (14c), is a precursor of the mentioned divergency. 

The structure factor or the pair correlation function de- 
termine the positions of the liquid-gas transition points. 
However, one faces the known consistency problem that 
different routes for the equation of state yield different 
results for the transition points if approximations for Sq 
or g{r) are used 0|. We will not discuss these problems 
in this paper, since it is irrelevant for the evolution of 
glassy dynamics or the glass transition whether the fluid 
is in a stable or metastable thermodynamic state. 



III. THE PHASE DIAGRAM 

A. The bifurcation equation 

The MCT equations of motion for various dynamical 
quantities are based on the equations for the normalized 
density correlators (l)q{t) = {p*^{i)Pq) / {\Pq\^) ■ For the fiq- 
uid state these functions approach zero for large time i; 
density fluctuations which were created at time t — Q 
disappear for t oo. The glass state is characterized by 
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a spontaneous arrest of these fluctuations, i.e. the long- 
time limits fq of the correlators do not vanish. The ideal 
liquid-glass transition of the MCT is characterized by a 
discontinuous increase of fq from its value zero in the 
liquid state to the critical Debye- Waller factor > of 
the glass. For colloidal suspensions, fq can be deduced 
from the dynamical-light-scattering results for (j>q (t) . The 
experimental findings for the HSS ||l0|| and for a charge 
stabilized system confirmed the discontinuity for fq 
and the data for agree well with the MCT results. 

The fq obey the equation - /,) = J^q{f) @. 

Here, the mode-coupling functional J-'q is given by 



'.^n^hff^f\n-k\ 



(17a) 



The mode-coupling vertices are determined by the struc- 
ture facto] 
density p 



ture factor 5*^, the direct correlation function Cg, and the 



\k = SqSkS^^y^^^P 



l<f-fcl 



(17b) 



In the following, the wave-vector integrals will be dis- 
cretized to points on a grid of M values, which are equally 
spaced with step size h, starting at gmin — h/2. Thereby 
the mode-coupling functional is changed to a second or- 
der polynomial 



^q{f) - X! ^Q^kpfkfp ■ 
kp 



(18) 



The explicit rera-esentation of the coefficients can be 
found in Ref. |^l|. The M parameters fq obey the al- 
gebraic equations 



fq/{l-fq)^Tqif). 



1, 



(19) 



Besides the long-time limit fq, Eq. (g^ can have fur- 
ther solutions, say fq, obeying < < 1. The Debye- 
Waller factor is distinguished by the maximum property 



fq > fq, q ^ 1, 



We used the iteration pro- 



cedure = .Fj/(")]/(l + .Fj/(")]), n ^ 0, 1, . . . 

to determine fq. With increasing n the decrease 
monotonically towards fq, if the iteration is started with 

/f = 1 El- 

Two concepts are needed in the following, namely 
the maximum eigenvalue E and the exponent parame- 
ter A [p2[. For the discussion of the implicit equations, 



Eq. (19), the Jacobian J is of importance. It is equiva- 
lent to 1 — C, where the M x M matrix C is determined 

by 



qk 



dTqjf) 

dfk 



a-fkf 



(20a) 



Also the variation of C with changes of / is needed 



Cq,kp — 



1 d^Tqif) 

2 dfkdfp 



{l-f,f il-f,f 



(20b) 



There is a nondegenerate eigenvalue E of matrix C with 
the property, that all other eigenvalues E obey E < E. 
There holds E < 1, and liquid-glass transition points are 
determined by the condition E — Ec — 1. It is helpful to 
follow the drift of E with changes of control parameters 
while searching for the transition points. The left and 
right eigenvectors of C for the eigenvalue E, denoted by e 
and e respectively, are uniquely determined by the condi- 
tions: eq > 0, e, > 0, Y.q ^qSq = 1, Eg 6,(1 - /g)e^ = 1. 

They are used to characterize every transition point by 
a single number A, defined as 



A 



E 

qkp 



'^q^ qjkp'^k^ 



(21) 



The solutions of Eq. (|l^), considered as functions of 
the coefficients Vq,kp, can exhibit singularities, which 
are called bifurcation points [Q . The singularities occur 
if the Jacobian J is a singular matrix, i.e. if the matrix C 
has eigenvalue unity. The special singularities, which are 
exhibited by the Debye- Waller factors, are called glass- 
transition singularities. These are members of the sim- 
plest family of singularities, labeled Ai, I = 2, Z, . . . [Q. 
They are topologically equivalent to the bifurcation sin- 
gularities of the real roots of real polynomials of degree 
I. Since the Vq^kp are smooth functions of the control pa- 
rameters, in the SWS the Debye- Waller factor fq exhibits 
Ai singularities considered as a function of the variables 
Lp, 9, and 5. The liquid-glass transition is an example 
for the simplest bifurcation singularity A2, called a fold 
bifurcation. Such transitions occur on smooth surfaces 
in the 3-dimensional parameter space. These surfaces 
can terminate in smooth lines of A3 singularities, that 
are also called cusp bifurcations. The inner points of 
the -surfaces are characterized by < A < 1, and 
for the endpoints there holds A = 1. The most compli- 
cated generic singularity in a three-parameter system is 
the meeting of two A^ lines in an Ai point. Its position 
shall be denoted by if* , 9* , and 5* . This singularity is 
also called a swallow-tail bifurcation . The possibility 
of the described scenarios has been demonstrated earlier 
for schematic MCT models, invented with the mere inten- 
tion of demonstrating the existence of A^ and A^ points 
This paper is the first demonstration of the exis- 
tence of an Ai for a microscopic model; the values for 
the SWS are given in Eq. (p^). 

The numerical work is done with step size hd = 0.4. 
It was checked for representative cases, that choosing 
smaller step sizes does not alter the results to be pre- 
sented. Choosing M is equivalent to introduc ing a n up- 



The 
were 



per wave- vector cutoff q* = {M — 1/2) /i in Eq. (|17a 
previous comprehensive studies for the HSS [ [liy45| 
done with M = 100. For sufficiently large 1/9 and suffi- 
ciently small S, the direct correlation function Cq develops 
a large-g tail, discussed in Eq. (I4h. This is decisive for 
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fixing the value q* needed to ensure th e cor rect handling 
of Eq. dl^) . If the coefficient K in Eq. ( |l4c|) is kept fixed, 
the cutoff q* will increase with decreasing well-width pa- 
rameter proportional to 1/S. The maximum value for M 
that can be handled in the numerical work, defines the 
lower limit for 1/K and 5, which can be treated. We used 
values for AI up to 2000 occasionally, in order to guar- 
antee the cutoff independence of the results reported in 
this paper. 



B. Results 

The phase diagram for the SWS is shown in Fig. |^ 
for several constant-(5 cuts through the three dimensional 
control-parameter space. The results based on the PYA 
and the MSA are qualitatively the same. Let us first 
consider the three states 1,2, and 3 from Fig. ^ for (5 = 
0.06. Within the MSA, state 1 refers to the liquid phase, 
cf. Fig. Increasing ip to the state 2 increases the 
height of the first sharp diffraction peak of Sq, located 
near qo = 7, cf. Fig. ^. Thus the compressibility for 
fluctuations in the shell q ^ qo increases. This leads to 
arrest in a glass state, as known from the HSS. If one 
cools state 1 at fixed Lp = 0.50 down to state 3, S, 
decreases, as was explained in connection with Fig 
This effect stabilizes the liquid, but it is overcompensated 
by the increase of Sq in the small-wave-vector region, 
q < 6.1, and, more important, in the large-g region, q > 
7.4. As a result of this compressibility increase on the 
wings of the structure factor peak, the liquid freezes to 
a glass upon cooling, cf. Fig. ^ For large temperature, 
Sq depends only weakly on T; the terms proportional to 
K ^ S ■ uo/iksT) in Eqs. (|) to (|ll|) cause only small 
modifications of the coefficients determining the factor 
function Q{r). This explains, why the transition lines 
are almost vertical in Fig. ^a for ksT/uQ > 3 and in 
Fig. I^b for 9 > 1. The peak wings are not very sensitive 
to density changes by a few percent for ip < 0.51. This 
explains, why the transition lines in Fig. ^ are rather flat 
there. The two pieces of the transition line join smoothly 
and A remains below unity for S = 0.06, as is shown 
in Fig. Thus, the described curve represents a cut 
through a smooth surface of A2 bifurcations. 

The mentioned high-temperature pieces of the tran- 
sition surface are located at packing fractions ipc, which 
exceed the value for the HSS, ip^^^ w 0.516. This means, 
that the attraction forces have stabilized the liquid phase. 
This effect is smaller for larger T and therefore the pc- 
versus- Tc curve decreases. There is the possibility of glass 
melting due to cooling, if the decrease of Sq„ is not over- 
compensated by the increase of the structure-factor-peak 
wings. The attraction causes bonding, in the sense that 
the average separation of two particles is smaller than ex- 
pected for a HSS. Therefore the average size of the holes 
increases and this favors the long-distance motion char- 
acteristic for a liquid. Consider states with p = 0.52 and 



5 = 0.06. For the MSA results one notices from Fig. ||b, 
that the system is in a glass for 6 — 0.10 and it melts 
upon heating if 9 approaches 9^^ w 0.30. This transition 
occurs due to the temperature drift of the coupling to 
modes with wave vectors in the wings of the Sq peak. 
The system remains in the liquid upon further heating 
until it reenters the glass at 6*+ w 1.41. This freezing is 
caused by the drift of the coupling to modes with q ^ qo. 
The described reentry phenomenon |^,^ is a manifes- 
tation of two mechanisms for localization due to the cage 
effect in the high density SWS to be explained below. 

The preceding two paragraphs can be summarized as 
follows. There is a subtle interplay of excluded- volume 
and bonding effects which determine the variations of the 
pair correlation function g(r) on the length scale of the 
particle diameter d. This is reflected in the properties 
of the structure factor Sq for wave numbers q at and 
around the peak position qo; the relevant q range is the 
one exhibited in Fig. ^j. Fluctuations with longer wave 
length are of no qualitative importance for those param- 
eter points studied in this paper. This conclusion was 
corroborated by dropping all contributions to the mode- 
coupling functional, Eq. ([T^, where k or p are smaller 
than 4; there was no significant change of the phase di- 
agram calculated with the MSA structure compared to 
what is shown in Fig. ||b. Similarly, a cutoff q^ = 20 
was introduced such that all contributions to the mode- 
coupling functional with q > q^ are dropped. The MSA 
results shown in Fig. for 5 — 0.09 and S = 0.06 did not 
change, nor was there a noticeable change for the other 
curves for 9 > 0.6. We conclude, that the two specified 
sources for correlations on intermediate length scale ex- 
plain the phase transition points, which are marked by 
open symbols in Fig. ^ or by the corresponding light 
lines in Fig. ||b. 

To substantiate the previous conclusion, we have con- 
structed a further phase diagram based on the MSA 
structure factor, using the above specified cutoff (74.. As 
mentioned, the curves for 6 — 0.09 and (5 = 0.06 were 
reproduced up to minor deviations. Upon decreasing b 
further, curves emerge, which continue the trend of the 
two ones for larger b. The limit (5 — > can be carried out; 
no new features appear. This ad-hoc MCT model yields 
a smooth liquid-glass transition surface of A2-bifurcation 
points with A < 1. Obviously, the true phase diagram 
of the SWS shown in Fig. || is quite different. There are 
the transition points marked by filled symbols in Fig. |^ 
or by the corresponding heavy lines in Fig. These 
define transition lines, which do not join smoothly the 
previously discussed lines. Rather they cross the former 
lines. Thus, for sufRciently small b and sufRciently large 
attraction strength uq/ (ksT), there is a new glass forma- 
tion mechanism, dominated by density fluctuations with 
large wave number g > g+ . These are due to spatial cor- 
relations on the length scale of the attraction-well width 
A. 

For a discussion of the identified new pieces of the 
transition surface, the mode-coupling coefficients in 
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Eqs. (17a) and ( |l7b| ) can be simplified. As explained 
in connection with Fig. |[ one can write for k = \k\ > q+ 
and p — \q — k\ > q+ the leading asymptotic expres- 



sion for the structure factors, 
direct correlations function c, 



Sk 



1, 



and for the 
Thus, 



, Eq. (|l4bD 

the dominant part of the mode-coupling functional de- 
pends explicitly on the control parameters by the pref- 
actor qK'^ = g {Suq/ (kBT))^ and otherwise only on 6 
via the large- wave- vector cutoff Qu- A further density 
dependence is due to fluctuations with q < only. 
This explains why the low-T-transition lines in Fig. 
are so flat. decreases proportional to 6'^, and this 
effect is not overcompensated by the increase of g„. As 
a result, the horizontal transition lines decrease with de- 
creasing 6. For low packing fractions, such a trend can 
be shown explicitely by an analytic calculation p5| , ^ . 
Along this transition line, A increases with increasing tp 
until it approaches unity signalizing an endpoint at some 
(y9°((5), 0°{S). This is demonstrated in Fig. ^ for three 
values of S and in Fig. ||b for S = 0.03. The line stops 
the previously discussed line in some crossing point. Be- 
tween the crossing point and the endpoint, there occur 
glass-glass transitions. Upon increasing S, the length of 
the glass-glass transition line shrinks, until it vanishes 
for some S* at some ip* — ip°{6*) and 9* — 9°{S*); and 
this is the A4 singularity whose coordinates are listed in 
Eq. @. 

Every pair of glass states can be connected by a curve 
in parameter space such that, upon shifting the control 
parameters (p, 9, and 6 along this curve, the glass proper- 
ties change smoothly. Thus one cannot discriminate pre- 
cisely between repulsion- and attraction-dominated glass 
states. However, upon crossing the glass-glass-transition 
surface, there occurs a discontinuous change of the glass- 
state properties. These discontinuities can be used to dif- 
ferentiate quantitatively between the two types of glasses. 
From a mathematical point of view, the situation is anal- 
ogous to the termination of the liquid-gas-transition line 
at the critical point. Some results shall be presented from 
the theory based on the PYA for Sq. Let us consider 
states for ip « 0.54 for the smallest S used in Fig. ||a. 
Figure exhibits as filled symbols the decrease of the 
Debye- Waller factor for three representative wave num- 
bers, if the attraction-dominated glass is heated towards 
the transition temperature Tc, ksTc/uQ — 1.0471 — 9c- 
The fq decrease towards the critical value according 
to the asymptotic square root law fq — fq o: ha\/9c — 9, 
which is the signature of the fold bifurcation E2]. Upon 
crossing the line into the repulsion-dominated glass, fq 
drops and keeps on decreasing upon further heating up 
to 9 Ki 1.3, as shown by the open symbols. Notice, that fq 
does not exhibit any singularity for 9 decreasing towards 
9c- The remarkable variation of fq for 9 near but above 
9c is a precursor phenomenon of the nearby ^3 singular- 
ity. The fq for 9 > 9c is smaller than the Debye- Waller 
factor for the HSS at the same packing fraction. Hence 
the attraction has softened the glass. This effect has to 



disappear for very large T, and this explains, why fq in- 
creases again, reflecting a glass stiffening upon heating. 
The described effects for 9 > 9c are the counterparts to 
what was discussed above in connection with the reentry 
phenomenon. 

The wave-vector dependence of fq changes qualita- 
tively upon crossing the glass-glass-transition line as is 
shown in Fig. ^ Let us focus in this paragraph on 
the wave- vector regime at and above the structure fac- 



tor peak position q 



> 



90 



7. Here, fq oscillates with 



wave- vector scale 2^/(1 around the Mofibauer-Lamb fac- 
tor /|. The latter is the analog to fq, constructed for a 
tagged particle with position vector (t) via its density- 
correlator </)^(i) = {pt{typt); pt{t) = exp (zgf;(0). This 
quantity, in particular its long-time limit can also 



be measured [ll3|. One finds f^ 



1 



'r^ + 0{q^), 



where r; is the localization length of the particle: — 
limt^oo(|^s(i) — 'rs(O)p) In the Gaussian approx- 

imation, one can write fq = exp {—q^rf) so that the 
half-width- wave vector qi, defined by /g, — 0.5, can be 
used to estimate r; « 1 /qi . For the Debye- Waller factors 
on the high-temperature side of the transition line, which 
are shown as full lines in Fig. ^, one estimates qi « 20. 
The localization length n = 0.05 is about the same size 
as expected for a particle rattling between the hard walls 
of its cage in a HSS. The fq and r; are close to those of 
a HSS [ p2p^ ]. However, on the low-temperature side of 
the transition, the fq vary by less than 10% if q increases 
up to 2qo, as shown by the dashed lines in Fig. ^. The 
wave number qi is much larger than expected from the 
free volume in the cage. At the transition, qi « 44, i.e. 
the localization length is decreased discontinuously by a 
factor of about 2.3. For 9 = 0.9 and 0.6 the localiza- 
tion length is about 0.01 and 0.006, respectively. This 
shows that the localization length is of order S- The par- 
ticle is bound to the wall of the cage and localization is 
determined entirely by the particle attraction. This lo- 
calization mechanism via bond formation is operative at 
low packing fractions also, and it has been studied within 
MCT in this regime. There, bond formation has been ar- 
gued to be of importance for colloidal gelation ||2^,^,^ . 

There is no strong dependence of fq on wave numbers 
g < 4. The low-temperature glass is distinguished from 
the high-temperature one by the fact that the zero- wave- 
number limit of the Debye- Waller factor, /g, is larger 
for the former than for the latter. Therefore, upon cross- 
ing the glass-glass-transition surface by cooling, the peak 
for (7 w go of the /g-versus-g diagram disappears. Fig. ||. 
The number /o is related to the longitudinal elastic mod- 
ulus of the system. This consists of a part expected for 
the ergodic liquid and a part toq reflecting the incom- 
plete relaxation of the non-ergodic glass |Q. The lat- 
ter is given by the zero-wave-number limit of the mode- 
couphng functional, Eq. (|l8|), mo = ^ojf) - The q = 
limit can be carri ed o ut easily in Ec[. (17b), so that one 
derives from Eq. (17a) a formula [n 
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Too = 



dkv^{k)fl 



(22) 



After discretization, one can substitute the results for /j. 
to get Too, and via Eq. ( p^ ) one has /o = Too / (1 + too). 
For the glass with 9 > 9c, the integral in Eq. ( P4 ) is dom- 
inated by fluctuations on the intermediate wave-vector 
scale, k « qq. As known from the HSS, toq is of order 
unity and thus /o is about 0.5. However, for 6 < 9c, 
the integral is dominated by large-fc fluctuations. This 
enhances too. For 9 = 9c, one finds too ~ 12.7 and this 
increase of the modulus explains the increase of /o to 
about 0.927, exhibited by the lowest dashed line in Fig. ||. 
Decreasing the temperature to ksT/uo = 0.6, leads to 
Too ~ 385 and this explains the large value /o ~ 0.997, 
exhibited by the uppermost dashed curve in Fig . ^. For 
the shear modulus G", a formula like Eq. (g2|) holds, 
where the expression for is similar to that for [ p^ . 
Therefore, the shear modulus also exhibits the specified 
strong enhancement due to the attraction wells. Figure ^ 
shows that the dramatic change of the moduli is the most 
relevant feature to be observed upon crossing the glass- 
glass-transition surface. The strong short-ranged attrac- 
tion causes bond formation, and this increases the rigid- 
ity of the glass with respect to compressions or shearing 
considerably relative to that of a glass at a similar den- 
sity, where the structural arrest is dominated by mere 
hard-sphere repulsion. 



IV. STRUCTURAL RELAXATION 

A. Some general MCT equations 

The MCT equations of motion are based on the exact 
expression of the density correlator (t)q{t) in terms of a 
fluctuating- force correlator jJJ^, which in turn is split 
into a part treating normal liquid effects, and a relaxation 
kernel mq{t) describing the cage-effect contribution. This 
kernel is approximated by the mode-coupling functional, 
discussed above in connection with Eqs. (17) and (18): 
j{t) = J'q [0(i)], [0j42|. This paper wiU be restricted to 



the simplest approximation for the normal liquid effects, 
i.e. the colloid will be treated as a system of Brownian 
particles, so that only the instantaneous correlations as 
given by the structure factor Sq are incorporated. As a 
result one obtains p| : 

TqdtMt) + Cl)q{t) + f mq(t-t')dt'<l>q{t')dt' ^0 . (23) 

Jo 

This equation implies the short-time asymptote (t>q{t) = 
1 — (t/Tq) + 0{t^). For the time scale, one finds Tq = 
Sql {D^q^), where denotes the single-particle diffusion 
coefficient. Dq reflects the property of the solvent, and 
it fixes the time scale for the transient motion. The unit 
of time shall be chosen such that l/Uo = 160 to ease 



comparisons with the results for the HSS from preceding 
work |l|,||. 

Two comments on the implications of Eq. (|2^) might 
be appropriate. First, the solutions are completely mono- 
tone functions, i.e. there is a rate density Pq{"f) > 0, nor- 
malized to pq (7) c?7 = 1 , such that 



(24) 



Thus, the MCT approximations maintain a fundamen- 
tal property of colloidal dynamics: auto-correlation 
functions can be written as superpositions of Debye- 
relaxation functions ]4^ . Second, outside the transient, 
the solutions can be written as <j)q{t) = Fq{t/to). Here, Fq 
is given by the mode-coupling functional Tq, i.e. by the 
equilibrium structure factor 5', . The transient dynamics, 
no matter how complicated, enters via the single time 
scale to only ||l6|,E8 49 1 . The following results for the long- 



time dynamics are thus not influenced by the simplified 
treatment of the short-time dynamics in Eq. (^3|) , except 
up to a change of the overall time scale tg. It is known 
that the short-time dynamics in colloids is influenced by 
hydrodynamic interactions |Q. Unfortunately, it is not 
known how to incorporate these interactions in a theory 
for high-density colloids. But we consider it plausible, 
that the hydrodynamic interactions merely renormalize 
the transient dynamics |5^], thereby being irrelevant for 
the structural-relaxation effects. 

For control parameters approaching a glass-transition 
singularity, there appears an increasingly larger dynam- 
ical window, where the solutions are arbitrarily close to 
the critical Debye-Waller factor Therefore, one can 
solve the MCT equations of motion by an asymptotic ex- 
pansion, using 5(j)q{t) = (jiqit) — fq as a small parameter. 
The result can be expressed in the form 



(25) 



Most of the known universal results for the MCT bifur- 
cation dynamics are based on the understanding of the 
leading-order contribution G{t). The next-to- leading- 
order term G*^^^ (t) allows to discuss the range of validity 
of the leading-order formulas. A comprehensive demon- 
stration of the cited results for the HSS can be found in 
Refs. ^ andpl. 

Equation (|25) reduces to a particular transparent form 
for the critical dynamics, i.e. for control parameters 
placed on the glass-transition points. For the fold bi- 
furcation, one finds a power-law decay: 



G(^) = (^oA)^ GW(t) = (ioA) 



2a 



A. 



(26a) 



The critical exponent o is given by the exponent parame- 
ter A via r(l - a) Vr(l - 2a) = A, < a < 0.5. The end- 
points of the 742-bifurcation surfaces are characterized by 
exponent a approaching zero. For the cusp bifurcation, 
there holds M: 
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G{t)^l/Mt/tof, 
G(i)(t)«ln(ln(tAo))/ln(tAo)' ; ^3 



(26b) 



For the swallow-tail bifurcation, one has plIJS 



G(t)(x l/ln(tAo), 
G'^^\t) (X ln(ln(i/io))/ln(i/io)^ ; A4 . (26c) 

The dependence of the leading-order contribution G{t) 
on the control parameters is well understood, but shall 
not be considered in this paper. 

B. Results 

Let us first estimate the dynamical window relevant for 
the discussion. For the HSS, the correlators (j>q{t) decay 
from 1.00 to 0.95 for times increasing up to about t = 1 
for representative wave vectors. In this sense, < = 1 is 
the scale for the transient dynamics. After a crossover 
window of about one or two decades, the leading-order 
asymptotic law 4>q{t) = fq + hqG(t) becomes valid at 
about t — 10^. This value may be an order of magnitude 
smaller or larger, depending on the wave number q [^. 
The same is true for the data obtained by van Megen 
et al. for hard-sphere colloids, provided one identifies 
the time unit t = 1 with 1 msec . The correlators 
have been measured up to 10® msec, and thus the so far 
explored windows extend up to 10^ in the units used here. 
This limit might shift up in future work, using different 
experimental setups. 

Figure |lO| exhibits the critical correlators 0^(i) for 
q = 4.2 for five states on the transition line 6 = 0.0465 
through the A4 singularity. The states 1 and 5 refer 
to an exponent parameter A ~ 0.80 on the side of the 
attraction-dominated and repulsion-dominated glass, re- 
spectively. For times of the order of 10^ and larger, the 
leading-order formula, (f>q{t) — oc {to/t)'^, a « 0.28, 
describes the results. Thus, the scenario is similar to 
the one known from the HSS, and this is also true for 
other states on the line with A < 0.80. However, if one 
considers states closer to the A4 point, the onset of the 
critical power law gets shifted to larger times. This is 
demonstrated for the two states 2 and 4, which deal with 
A « 0.90. For the state 4, the t'" law with a « 0.20 
is valid only for t > 10®. This trend continues, if one 
moves even closer to the A^ point, whereby A increases 
even further; compare Fig. ^. At the A4 singularity, the 
correlator decays from 1.00 to fg ~ 0.77. This decay is 
stretched so enormously, that even for t — 10^^ it only 
reaches the value 4>q{t) ~ 0.80, as is shown by curve 3 
On e ca n describe the critical correlator with 
and (p6^) for the window 10^^ < t < 10^^, us- 




ing hq and hq 



(1) 



as fit parameters. But this fit does not 



describe the correlator for t < 10^". Thus one concludes, 
that the critical correlator of the A4 for the SWS cannot 



be described by the asymptotic Eq. ( 26c ) within accessi- 
ble dynamical windows. Nor can the critical power-law 
decay of the A2 singularity be measured, if A exceeds a 
certain value, say 0.9. Thus, there is a part of the tran- 
sition lines near the A4 point, characterized by A > 0.9, 
where the correlators exhibit structural relaxation pat- 
terns towards the plateau values /| that are stretched 
up to t = 10® or larger. The known asymptotic formulas 
cannot be used to describe the MCT solutions within this 
regime. 

The liquid dynamics on the small-(5 side of the A4 point 
is particularly subtle, since there is an A3 singularity in 
addition to the line-crossing point. Figure |ll| exhibits as 
an example such a situation for 6 = 0.03. Parameters on 
a straight line, 9 = 0.1875, which is slightly above the 
A3 point, are considered. The transition to a repulsion- 
dominated glass state then occurs at ipc — 0.5360. At 
the transition point, the critical Debye- Waller factor is 

(l]c 

fq ~ 0.50, and the exponent parameter is given by 
A = 0.847, implying a critical exponent a = 0.250. Curve 
3 was calculated for such small distance from the tran- 
sition point, — £ = {ific — ^)/^c = 7.9 • 10""', that 0g(i) 
decays to zero just within the dynamical window dis- 
played in the figure. The dash-dotted line with label A 
presents the leading-order critical law for the A2 singu- 
larity, Eqs. ( p5| ) and ( ^6a| ) , with the time scale to fit- 
ted to the decay at long times for ip = (pc- One ob- 
serves the same phenomenon as explained above in con- 
nection with Fig. |l^: since A is rather large, the asymp- 
totic law describes the data only for rather large times, 
t > lO^ '"^. After falling below the plateau value fq^^", 
the correlator decays towards zero. This is the a pro- 
cess, and it starts with the von-Schweidler asymptote 
fq^^"^ — hqit/TY , which is shown by the dash-dotted hne 
with label B. The exponent h = 0.396 obeys a similar 
relation as the critical exponent, r(l-|-&)^/r(l-|-26) = A. 
Thus the structural relaxation connected with the liquid- 
glass transition follows the known scenario, except that 
the familiar A2 patterns can be observed only for times 
far out the transient regime, t > 10^'^. 

Figure |l^ exhibits a large dynamical window, 10^ < 
t < 10^-^, where the structural relaxation does not fol- 
low the asymptotic laws for a fold bifurcation. Instead, 
the dashed straight line demonstrates that the correlator 
labeled 3 follows a logarithmic decay law, 



jit) 



/•(2)c 
J a 



Gqlnt 



(27) 



for the major part of the mentioned window, 10^ < t < 
lO'*-^. Here, /f « 0.87 is the Debye-Waller factor at 
the A3 singularity. There is a line through every A3 sin- 
gularity, which is transversal to the transition line ending 
at the A3, such that Eq. ( |2^ ) is a leading-order solution 
for the MCT equations of motion on a certain interme- 
diate time window. The length of the window expands 
and the prefactor Cq in front of the Ini decreases, if one 
moves closer towards the A3 point Q . These results ex- 
plain the appearance of the In t part and the change of its 
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prefactor, if one compares curve 2 with curve 3 in Fig. |T^. 
One concludes that it is the bifurcation dynamics of the 
A3 singularity which prevents the evolution of the 
law for the fold bifurcation. Similarly, the a process for 
curve 2 does not start with von Schweidler's law. There- 
fore, contrary to what one observes for the dynamics of 
the HSS for comparable large times the a processes 
for curves 2 and 3 do not obey the superposition prin- 
ciple. Close to the A3 point, the dynamics outside the 
transient and preceding the onset of the Int-decay law 
follo ws th e critical law for the A3, as given by Eqs. ( |25| ) 
and (26b). But for the shown curves, the situation is 
similar as explained comprehensively for the critical de- 
cay for the HSS |4l]. The plateau fjf'''^ is so high, that 
there is only a small variation remaining for the l/(lnt)^ 
law to manifest itself. The correction terms for the cited 
leading-order and next-to-leading-order contributions are 
so large, that one has to consi der s tates much closer to 
the A3 to see the result of Eq. ( 26b ). 

For state 1, the A2, A3, and A4 singularities are so far 
away, that none of the cited asymptotic laws is clearly 
developed. On the other hand, they are close enough to 
cause a considerable relaxation stretching. The correla- 
tor 4>q(t) needs a dynamical window of three orders of 
magnitude to complete the 80% of its decay from 0.9 to 
0.1, as is shown in Fig. |ll|. 

In Fig. [l^a, a set of correlators (j)q{t) for five repre- 
sentative wave numbers is shown. The state refers to 
the liquid close to the A4 point. For q = 24.2, the 
window for the logarithmic decay extends from t = 10'^ 
to t = 10^°. For the other wave numbers, the corre- 
sponding windows are smaller, i.e. the window of valid- 
ity for the leading-order asymptotic laws depends on the 
chosen correlator. If the correction term hq^'' G^^'^ (t) in 
Eq. ( ^5|) could be neglected, i.e. if the factorization theo- 
rem (t)q{t) — fq = hqG{t) would hold, the rescaled correla- 
tors (jiqit) = [(pqit) — ff\ /hq should coUapsc on the com- 
mon function G{t). In particular, all correlators should 
cross their plateau value at the same time given by 
G(t_) = 0. The latter property is fulfilled within a small 
error margin for t_ — 8.9 x 10^. Figure [l^ demonstrates 
the validity of the factorization property for a two-decade 
window. The size of this window is considerably smaller 
than the one found for the HSS for a state with a similar 
overall relaxation time El. Thus, the next-to-leading- 
order correction in Eq. p5[ ) is much larger near the ^4 
than the one known from the bifurcation dynamics of the 
simple HSS. 



V. CONCLUSIONS 

In this paper, ideal liquid-glass transitions and the evo- 
lution of glassy dynamics were analyzed within the basic 
version of the mode-coupling theory (MCT) for a simple 
colloid model, where the particles interact via a square- 
well potential. The discussion was restricted to the high- 



density regime. Hence the excluded-volume effects play 
a crucial role for the structure, and the cage effect is 
an essential feature of the dynamics. The presence of 
short-ranged attractions leads to a variety of new fea- 
tures compared to the ones known from the hard-sphere 
system (HSS). We find a subtle phase diagram for the 
glass-transition lines in the plane spanned by the two 
control parameters, packing fraction (p and reduced tem- 
perature 6 (Fig. H). The diagram is organized around 
an A4-glass-transition singularity. This occurs for a crit- 
ical value 5* k, 0.04 of the ratio 5 of the attraction-well 
width and the hard-core diameter, a packing fraction Lp* 
exceeding the transition density f^^^ of the HSS, and a 
certain critical temperature, cf. Eq. (p^. 

For 5 > 5* and sufficiently low temperature, there is 
a part of the liquid-glass transition line, where the criti- 
cal temperature 9c increases with the critical density 
As expected for conventional liquids with, e.g., Lennard- 
Jones interactions, the glass transition can occur either 
upon cooling or upon compression. This part of the tran- 
sition line extends up to densities where ip exceeds t/?^^^, 
since the bonding effects due to the attraction stabilize 
the liquid phase. For large temperatures, the effects of 
the attraction get suppressed. Therefore, there exists a 
high-temperature piece of the transition line, where 9c 
decreases with increasing pc- There appears a regime of 
high density, where the liquid can transform to a glass 
either by cooling or by heating. The possibility of such a 
reentry phenomenon is characteristic for systems with a 
hard-core repulsion. In a conventional system, the effect 
cannot occur, since a soft-core repulsion implies a de- 
crease of the effective repulsion-core diameter with heat- 
ing; and this decrease overcompensates the effect of the 
decrease of the effective attraction strength. 

For 5 < 5* , the two mentioned transition- line parts no 
longer join smoothly. Rather the low-temperature line 
terminates the high-temperature one at some crossing- 
point, such that they appear as two separate transi- 
tion lines. At very high temperature, the mechanism 
of glass formation is similar to the one of the HSS, 
and in general, the temperature dependence of the high- 
temperature transition line is weak. Glass transitions 
across this line are caused by an arrest of density fluctu- 
ations on the length scale of the inter-particle distance. A 
tagged particle is localized due to repulsion by its cage- 
forming neighbours. In contrast, the low-temperature 
line describes glass formation due to the arrest of den- 
sity fluctuations on a length scale of the order of the 
attraction-shell width. Here, tagged particles are local- 
ized due to a formation of short bonds with their cage- 
forming neighbours. The density dependence of these 
transition points is weak, and the transition line extends 
into the regime of gel formation at low densities. On the 
high-density side, it extends into the glass regime, until 
it ends at an Aa-glass-transition singularity, as indicated 
by the open circles in Fig. ^d. 

Between the mentioned line-crossing point and the 
endpoint of the second transition line, there is a line of 
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glass-glass transitions. The averaged equilibrium struc- 
ture, as characterized by the structure factor Sq^ is the 
same on either side of this line. But the two differ- 
ent localization mechanisms imply qualitatively different 
frozen structures, reflected by differences in the Debye- 
Waller factor fq. The one on the high-temperature side, 
shown by the uppermost solid line in Fig. ^, is similar to 
the Debye- Waller factor of the HSS at the same density. 
It exhibits a pronounced peak near the position go of the 
structure-factor peak, and the zero-wave number limit /q 
is about 0.7. On the low-temperature side, is much 
larger, as is shown by the lowest of the dashed lines in 
Fig. ^. In particular, /g is considerably increased. The 
/^-versus-g curve for the attraction-dominated glass is 
bell-shaped like a M66bauer-Lamb factor. The increase 
of /o towards the upper limit unity is connected with a 
large enhancement of the longitudinal modulus. Cross- 
ing the glass-glass transition line, the longitudinal mod- 
ulus as well as the shear modulus experiences a large 
discontinuity, as shown in Fig. ^. The large differences in 
the macroscopic elastic properties are the most obvious 
manifestations of the two localization mechanisms in the 
high-density system predicted by our theory. 

Two general MCT predictions for the relaxation near 
a critical temperature or critical density have been con- 
firmed by many experiments and molecular-dynamics 
simulations | |5^ . First, the structural relaxation exhibits 
a two-step scenario. Outside the transient, there oc- 
curs a relaxation towards the plateau For this step, 
cP(j)q{t) /d{lnt)'^ is positive. Then there is the a process 
dealing with the relaxation from the plateau to zero. Its 
initial part exhibits a negative second derivative of the 

(i)-versus-ln i curve. Second, there holds the super- 
position principle for the a process. On a time window 
that expands with increasing relaxation time, the 4>q{t)- 
versus-lnt curves can be collapsed on a common master 
curve by shifts along the abscissa. These two simple re- 
sults, which are fingerprints of the A2 bifurcation, are not 
valid for the relaxation at states close to an A4 singular- 
ity. The curves in Fig. ^ cannot be rescaled onto an 
a-relaxation master curve. The results in Fig. 03 do not 
exhibit changes of the second derivatives for ^^t) near 
the plateau fq. It was shown that the higher-order glass- 
transition singularities A3 and A4 cause strong pertur- 
bations of the asymptotic laws usually considered, valid 
close to the A2 bifurcation. In the present case, they 
can only be observed in windows, that might be outside 
the regimes accessible by experiments. In addition, the 
known asymptotic laws for the relaxation near A3 or A4 
glass-transition singularities also show up only in win- 
dows, that are irrelevant for experimental studies. These 
predictions of our theory do not seem to be a peculiarity 
of the square-well system. Similar results already hold 
for simple one-component schematic models 

An exception to the findings summarized in the pre- 
ceding paragraph is the logarithmic-decay law, Eq. ([27[). 
This characteristic feature of the dynamics near higher- 



order glass-transition singularities could be identified eas- 
ily in our results, as shown in Fig. Indeed, it is shown 
in Fig. |l^, that this In t decay is a precursor phenomenon, 
hindering the evolution of the A2 asymptotics. In par- 
ticular, there can be a crossover from the hit decay to 
the von-Schweidler decay around the point, where the 
(/ig (i)-versus-ln t curve crosses the plateau as is shown 
by curve 3 in Fig. |ll]. A similar scenario was recently 
observed for relaxation in a micellar system |55[ . 

The found extreme stretching phenomena have impor- 
tant implications for the experimental tests of MCT. In 
an experiment, it is not easily possible to measure self- 
averaged correlation functions for states like the ones dis- 
cussed in Fig. |l^. Thus, experimental results are likely 
to refer to history-dependent non-equilibrium states, and 
ageing effects are likely to be more pronounced than they 
are for the normal liquid-glass transition. Even if proper 
averaging could be achieved, one cannot determine the 
Debye- Waller factor fq — (j)q{t 00) within accessible 
time windows, if the states are similar to the ones shown 
with labels 3 and 4 in Fig. |o|. Similar conclusions apply 
for the measurements of the moduli near the glass-glass- 
transition line. 

The presented theory is based on some assumptions 
which we would like to discuss in more detail. First, 
one should expect that the equilibrium state of the sys- 
tem in the density regime considered is a crystal rather 
than the assumed amorphous phase. In experiments for 
colloids, crystallization is bypassed by chosing a poly- 
dispersity p for the particle diameters. Since nucleation 
rates decrease dramatically with increasing p, a choice of 
p of some percent is sufficient to establish a metastable 
amorphous state for practically arbitrarily long times. A 
small p causes only small changes of the calculated struc- 
ture factors, and thus only small changes in the coupling 
coefficients entering the MCT equations. Hence a small 
p will only imply small changes of the presented results. 
Indeed, it was shown for the HSS, that a change of p 
did not yield detectable changes of the measured 0g(t) 
p6| , ^ . But, it is unclear how strongly, e.g., the calcu- 
lated value 6* for the attraction-well width at the A4 
singularity will change, if a realistic value for p is consid- 
ered. 

The structure factor Sq of the stable or metastable 
equilibrium is used as input information for our work. 
Thus, the second source of reservations is due to the er- 
rors hidden in the used Sq. A well-known problem is 
that of the so-called thermodynamic inconsistency. Ther- 
modynamic quantities calculated along different routes 
using an approximate Sq as input often are not consis- 
tent with each other. Sophisticated closures involving 
adjustable parameters could be used to overcome this 
problem [Q. Alas, since thermodynamics deals with the 
q limit, for which the phase volume in the mode- 
coupling integrals becomes small, one would gain no fur- 
ther insight carrying out our calculations of Sees. Ill and 



[V using an improved closure for Sq. For the HSS, one 



finds only minor changes in the numerical values for the 
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transition points | |5^ , |5q | , and the same is anticipated for 
the SWS. A further difficulty arises regarding the smaU-r 
hmit of c(r) and g{r). Due to the approximations intro- 
duced for Q{r), one cannot guarantee that the excluded- 
volume effect, g{r < d) — 0, is exactly reproduced. In 
fact, we find that c(r) and thus g{r) develop a pole A/r, 
cf. Eq. (|l2|). Since g{r) is a distribution and since A/r 
is integrable in three dimensions, an A/r-teim is to be 
viewed as small, provided A is small. In the original work 
on the sticky hard spheres ||2^, A = 0{K'^(p). In our so- 
lution, A = 0{KSip). The hmits r ^ and 5 ^ do not 
commute, and our analytical solution decreases the error 
from a (5° to a 6^ effect. Since our results based on the 
Percus-Yevick closure and on the mean-spherical approxi- 
mation are in semi-quantitative agreement, we anticipate 
that better theories for Sq will not change the qualitative 
results of our theory. 

Third, the range of applicability of the MCT is not un- 
derstood. One can use the successful tests of the theory 
by the experiments performed on hard-sphere colloids, 
that were cited in Sec. |, as an a posteriori justification 
of MCT. But it is not clear, whether or not this theory 
can handle the effects caused by the formation of strong 
short bonds. On the other hand, the phenomenon of liq- 
uid stabilization due to bond formation and the resulting 
reentry effect, as well as the drastic changes of the elastic 
properties at the glass-glass transition, seem very plausi- 
ble. The fact that MCT brings out these subtleties might 
be considered as an argument in favour of this approach. 
In summary, it is the intention of this paper to point out 
the possibility of new features of glassy dynamics and 
to suggest a search for these features by experiments on 
colloids. 
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APPENDIX A: THE MSA FACTOR FUNCTION 

1. General Formulae 

Starting from Eqs. (^ and (B), Eq. for the factor 
function, and expressions Eq. ( ^0| ) and Eq. (^ij) for the 
parameters a, b, c shall be derived. The region < ?' < 
1 + S, for which Q{r) is nonzero, can be split into three 
parts: 



qi{r), <r < S 
Q{r) = { giiW, 5 <r <l 

qiii{r'), 0<r' = r-\<5. 



(Al) 



Q{r) is continuous at the boundaries of the intervals, in 
particular Q{1 + 5) = 0. From Eq. (||) together with 
g{r) = \+h{r) and G{r') = {l+r')g{l+r') the derivatives 
for the three parts of the factor function are obtained 



q[{r) = ar + b — 12(^ / ds G{s — r)qiu{s) , 
q'liir) =ar + b, 



(A2a) 
(A2b) 



q'uiir') = ar + b- G{r') + 12(^ / ds G{r' - s)qi{s) . (A2c) 

Jo 

Here g{r) = for < r < 1 was used, and the definition 
for a and b is given in Eq. (|6b|). The integrated form 
of En. (^ is used to introduce the MSA closure as in 
Eq. I 

qiii{r') = 12lp / dsqiii[s)qi{s - r') 



-K 



r' 5 f r 



(A3) 



where K — uoS/ksT. In the following, Eqs. (A2) and 
( |A3| ) are solved together with the reformulated expres- 
sions for a and b 



a = l-l2ip 
b = 12(p 



6 pi pS 

dsqi{s)+ / dsqu{s)+ / ds qiu{s) 
J 5 Jo 



5 /• 1 p5 

dssqi{s)+ / dssqii{s)+ / ds s qiii{s) 
J 5 Jo 



dsqni{s) 



Equation (A2b) gives Eq. (9a), 

<Zii(^) = ar^/2 + br + c 



(A5) 



where the continuity qii{r = 1) = = 0) yields c. 

For intervals I and III an expansion in d for fixed K will 
be performed. 



2. Leading Order 

In Eq. ( |A3|) , r' /S is of order (5° and the integral is of 
higher order, S^. Therefore, in leading order. 



qui{r') = K{l-r'/S) 



(A6) 



The boundary condition mentioned above fixes c== K — 
a/2 — b in leading order. Substituting Eq. (A6) into 
Eq. ( A2cD and keeping only terms in lowest or der i n S 
results in G(r') = K/S. The integral in Eq. (A2a) is 



(A4a) 



(A4b) 
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again of higher order and the continuity at the boundary 
r = S gives the parabola also for region I, 



1 ' r 



gi(r) ~ + br + c 



(A7) 



Inserting the factor function into Eq. (A4) and keeping 
only lowest order terms, linear equations for the param- 
eters are obtained 



a = 1 - 12(^(a/6 + 6/2 + c) 
6= l2ip{a/8 + b/3 + c/2) , 



which leads to Eq. (nU 



(A8a) 
(A8b) 



3. Next-To-Leading Order 

Substituting the leading order results into Eqs. ( [A^ ) 
and ( A3) pro duces the next- to-leading order. Subtract- 
ing Eq. (A2b) from Eq. ( A2a), the result for the interval I 
is given as the derivative 



q'j{r) - q'jj{r) = -12tp / dsK/6K{l - s/S) = 

J r 

= -6ipK^ {l-r/df , (A9) 

whic h is integrated to give the last term in Eq. (^). In 
Eq. (A3), the entire last line has to be taken into account 
for the next-to-leading order. The integral reads 



0.5 



0.4 



0.3 



0.2 



1 O ♦ 2 



40 ♦S 




9 

FIG. 1. Control-parameter plane for the square-well 
system (SWS), plotted as dimensionless temperature 
= kBT/uo versus packing fraction tp. The full line shows 
the spinodal calculated within the MSA for the relative at- 
traction-well width 5 — 0.05. Dashed (dash-dotted) curves 
show the corresponding spinodals for S = 0.03 {5 = 0.09). 
Vertical dotted lines mark the region for which the phase di- 
agram is discussed below in Fig. H. Diamonds mark the state 
parameters for which the structure factor is shown in Fig. 



12.^/ dsqiu{s)qi{s ~ r') ^ 

6K5-(pco{l-r'/Sf +0{6^) 



(AlO) 



Combinig Eqs. (A3) and ( |AlCl| ) yields the expression for 
the next-to-leading order term for the factor function in 
the outer shell, Eq. (|9^. The continuity at ?' = 1 in- 
troduces a modification of c from its leading order value 
Co, 



c^K-a/2-b + S- K/2 + 6(5 • Ktpco 



(All) 



where a and b are given by inserting the factor functions 
into Eq. ( |A4|) , 



a = l-l2ip (a/6 + b/2 + c + 5- K/2) 
h ^ 12ip (a/8 + 6/3 + C/2 + 5- K/2) . 



(Al2a) 
(A12b) 



This yields Eqs. 



(hop, and (M). 
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FIG. 2. Structure factor Sq and pair-correlation function 
g{r) of the SWS calculated within the MSA for relative well 
width S — 0.05. The labels 1 to 4 correspond to the states 
indicated by the diamonds in Fig. ^ They are given by the 
pairs {if, 6) of packing fraction and reduced temperature (0.50, 
0.50), (0.55, 0.50), (0.50, 0.15), and (0.40, 0.15), respectively. 
Here and in the following figures, the hard core diameter is 
chosen as the unit of length, d = 1, and q is given here and 
in the following in units of dT^ . 




FIG. 3. Direct correlation function Cq of the SWS for rel- 
ative well width 5 = 0.05, calculated within the MSA (solid 
line). Density and temperature are the ones considered in 
Figs. and |^ for the label 3. The dashed line exhibits the 
leading asymptote according to Eq. (|l^). The dotted line 
represents the same result with coefficient C replaced by zero 
(see text). 




-1.5 



FIG. 4. Factor function Qir) of the SWS. The hard-core 
diameter is chosen as the unit of length. The dashed line is 
the PYA result for 93 = 0.5293, reduced temperature 6 = 1.10, 
and well width parameter 5 = 0.0429. The full line is the 
MSA result for ip = 0.5258, 6 = 0.2332, 5 = 0.0465, chosen to 
represent the same physical state of interest in our discussion; 
see text for details. The dotted line shows the result for the 
HSS at packing fraction ip = 0.516. 
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FIG. 5. The phase diagram of the SWS showing cuts 
through the control parameter space for fixed relative attrac- 
tion-well width S = A/d. The upper part (a) is based on 
the PYA for the structure factor Sq, and the ratio (5/(1 + S) 
is noted in the legend. The lower part (b) is based on the 
MSA for Sq, and the weU widths are S = 0.09, 0.06, 0.0465, 
0.035, and 0.03, subsequently. The A3 endpoints are marked 
by open circles and the A4, by an asterisk. The vertical dashed 
line marks the transition line ip « 0.516 for the hard-sphere 
system. For reference, states 1 and 2 from Fig. |l|are included 
as diamonds. 
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FIG. 6. Exponent parameter A for points on three transi- 
tion lines. The upper part (a) was calculated within the PYA 
for the ratios S/{1 6) noted in the legend. Part (b) shows 
the results for the MSA, where the dashed line indicates the 
value A = 0.735 of the HSS. 



17 




FIG. 7. The Debye- Waller factor /, as a function of re- 
duced temperature for fixed packing fraction ip — 0.539672 
and fixed (5/(1 + 5) = 0.03. The wave vector qd = 7.0 is 
close to the structure factor peak position. The calculations 
are based on the PYA for Sq. The path through the pa- 
rameter space deals with a glass-glass transition occuring at 
8c = 1.0471, compare Fig. ||a. 




FIG. 8. The Debye-Waller factors for ip = 0.539672, 
5/{l + S) — 0.03. The temperature increases from top to 
bottom as 6* = ksT/uo = 0.600, 0.900, 1.000, 1.035, 1.0471 
(dashed lines) and 9 = 1.047, 1.150, 1.300 (full lines). The 
calculations are based on the PYA for Sq. 




FIG. 9. The shear modulus G' and the dimensionless lon- 
gitudinal elastic modulus mo as a function of the reduced 
temperature. The result is based on the PYA for Sq, and 
the parameters of state deal with the same path through the 
glass-glass transition as discussed in Figs. M and H. 
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FIG. 10. The full lines are the critical correlators (j>q{t) for 
the wave vector q — 4.2 calculated with the MSA-structure 
factor for the critical attraction well width S* — 0.0465. The 
states are located on the transition line as shown in the inset 
and refer to critical Debye-Waller factors fq = 0.973, 0.910, 
0.772, 0.593, 0.460 (from top to bottom as indicated by hor- 
izontal straight lines). For states 1 and 5, A = 0.80 corre- 
sponding to a critical exponent a = 0.279; state 2 (4) refers 
to A = 0.895 (0.908) corresponding to a = 0.210 (0.202). The 
dash-dotted lines show the asympto tes -\- Aqt^"" . State 3 is 
at the Ai singularity, given by Eq. (|l6b|) . 
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FIG. 11. Correlators ^q{t) for q = 4.2 calculated for the 
MSA-structure factor of a SWS with attraction-weU width 
5 = 0.03 and the reduced temperature 9 = 0.1875 for the three 
packing fractions ifi = 0.5000, ip2 = 0.5300, (fia = 0.5357 (full 
lines). The inset shows the relevant section from the phase 
diagram of Fig. |(b). The dashed-dotted lines with labels 
A and B exhibit the critical law /^^^ " + Ag/t°-'^^° and the 
von-Schweidler law fq^"^ — Bgt°'^^^ , respectively. The straight 
dashed lines exhibit logarithmic decay laws, Eq. (UTI), see text. 




FIG. 12. Correlators for a liquid state close to the 
A4 singularity calculated with the MSA for the struc- 
ture factor. The curves in the upper panel (a) show the 
(ttqit) where the labels 1 to 5 indicate the wave numbers 
q = 4.2, 7.0, 8.2, 20.2, 24.2. The corresponding critical De- 
bye- Waller factors flj are 0.764, 0.943, 0.860, 0.507, 0.369, re- 
spectively. The curves in the lower panel (b) exhibit the 
rescaled correlators (f>q{t) = [<^q(t) — /q] /hq. Here, the criti- 
cal amplitudes hq — {i — fq)^eq have the values 0.4665, 0.1343, 
0.2881, 0.7291, 0.7835. The dashed vertical line marks the 
time t- = 8.9 x 10^. 



